SETD3 regulates endoderm differentiation of mouse embryonic stem cells through canonical Wnt signaling pathway.

With self-renewal and pluripotency features, embryonic stem cells (ESCs) provide an invaluable tool to investigate early cell fate decisions. Pluripotency exit and lineage commitment depend on precise regulation of gene expression that requires coordination between transcription (TF) and chromatin factors in response to various signaling pathways. SET domain-containing 3 (SETD3) is a methyltransferase that can modify histones in the nucleus and actin in the cytoplasm. Through an shRNA screen, we previously identified SETD3 as an important factor in the meso/endodermal lineage commitment of mouse ESCs (mESC). In this study, we identified SETD3-dependent transcriptomic changes during endoderm differentiation of mESCs using time-course RNA-seq analysis. We found that SETD3 is involved in the timely activation of the endoderm-related gene network. The canonical Wnt signaling pathway was one of the markedly altered signaling pathways in the absence of SETD3. The assessment of Wnt transcriptional activity revealed a significant reduction in Setd3-deleted (setd3∆) mESCs coincident with a decrease in the nuclear pool of the key TF ß-catenin level, though no change was observed in its mRNA or total protein level. Furthermore, a proximity ligation assay (PLA) found an interaction between SETD3 and ß-catenin. We were able to rescue the differentiation defect by stably re-expressing SETD3 or activating the canonical Wnt signaling pathway by changing mESC culture conditions. Our results suggest that alterations in the canonical Wnt pathway activity and subcellular localization of ß-catenin might contribute to the endoderm differentiation defect of setd3∆ mESCs.

Epigenetic-focused CRISPR/Cas9 screen identifies (absent, small, or homeotic)2-like protein (ASH2L) as a regulator of glioblastoma cell survival.

BACKGROUND: Glioblastoma is the most common and aggressive primary brain tumor with extremely poor prognosis, highlighting an urgent need for developing novel treatment options. Identifying epigenetic vulnerabilities of cancer cells can provide excellent therapeutic intervention points for various types of cancers. METHOD: In this study, we investigated epigenetic regulators of glioblastoma cell survival through CRISPR/Cas9 based genetic ablation screens using a customized sgRNA library EpiDoKOL, which targets critical functional domains of chromatin modifiers. RESULTS: Screens conducted in multiple cell lines revealed ASH2L, a histone lysine methyltransferase complex subunit, as a major regulator of glioblastoma cell viability. ASH2L depletion led to cell cycle arrest and apoptosis. RNA sequencing and greenCUT&RUN together identified a set of cell cycle regulatory genes, such as TRA2B, BARD1, KIF20B, ARID4A and SMARCC1 that were downregulated upon ASH2L depletion. Mass spectrometry analysis revealed the interaction partners of ASH2L in glioblastoma cell lines as SET1/MLL family members including SETD1A, SETD1B, MLL1 and MLL2. We further showed that glioblastoma cells had a differential dependency on expression of SET1/MLL family members for survival. The growth of ASH2L-depleted glioblastoma cells was markedly slower than controls in orthotopic in vivo models. TCGA analysis showed high ASH2L expression in glioblastoma compared to low grade gliomas and immunohistochemical analysis revealed significant ASH2L expression in glioblastoma tissues, attesting to its clinical relevance. Therefore, high throughput, robust and affordable screens with focused libraries, such as EpiDoKOL, holds great promise to enable rapid discovery of novel epigenetic regulators of cancer cell survival, such as ASH2L. CONCLUSION: Together, we suggest that targeting ASH2L could serve as a new therapeutic opportunity for glioblastoma. Video Abstract.

Neurodevelopmental disorders and cancer networks share pathways, but differ in mechanisms, signaling strength, and outcome.

Epidemiological studies suggest that individuals with neurodevelopmental disorders (NDDs) are more prone to develop certain types of cancer. Notably, however, the case statistics can be impacted by late discovery of cancer in individuals afflicted with NDDs, such as intellectual disorders, autism, and schizophrenia, which may bias the numbers. As to NDD-associated mutations, in most cases, they are germline while cancer mutations are sporadic, emerging during life. However, somatic mosaicism can spur NDDs, and cancer-related mutations can be germline. NDDs and cancer share proteins, pathways, and mutations. Here we ask (i) exactly which features they share, and (ii) how, despite their commonalities, they differ in clinical outcomes. To tackle these questions, we employed a statistical framework followed by network analysis. Our thorough exploration of the mutations, reconstructed disease-specific networks, pathways, and transcriptome levels and profiles of autism spectrum disorder (ASD) and cancers, point to signaling strength as the key factor: strong signaling promotes cell proliferation in cancer, and weaker (moderate) signaling impacts differentiation in ASD. Thus, we suggest that signaling strength, not activating mutations, can decide clinical outcome.

Neurodevelopmental disorders, like cancer, are connected to impaired chromatin remodelers, PI3K/mTOR, and PAK1-regulated MAPK.

Neurodevelopmental disorders (NDDs) and cancer share proteins, pathways, and mutations. Their clinical symptoms are different. However, individuals with NDDs have higher probabilities of eventually developing cancer. Here, we review the literature and ask how the shared features can lead to different medical conditions and why having an NDD first can increase the chances of malignancy. To explore these vital questions, we focus on dysregulated PI3K/mTOR, a major brain cell growth pathway in differentiation, and MAPK, a critical pathway in proliferation, a hallmark of cancer. Differentiation is governed by chromatin organization, making aberrant chromatin remodelers highly likely agents in NDDs. Dysregulated chromatin organization and accessibility influence the lineage of specific cell brain types at specific embryonic development stages. PAK1, with pivotal roles in brain development and in cancer, also regulates MAPK. We review, clarify, and connect dysregulated pathways with dysregulated proliferation and differentiation in cancer and NDDs and highlight PAK1 role in brain development and MAPK regulation. Exactly how PAK1 activation controls brain development, and why specific chromatin remodeler components, e.g., BAF170 encoded by SMARCC2 in autism, await clarification.

Pan-cancer clinical impact of latent drivers from double mutations.

Here, we discover potential 'latent driver' mutations in cancer genomes. Latent drivers have low frequencies and minor observable translational potential. As such, to date they have escaped identification. Their discovery is important, since when paired in cis, latent driver mutations can drive cancer. Our comprehensive statistical analysis of the pan-cancer mutation profiles of ~60,000 tumor sequences from the TCGA and AACR-GENIE cohorts identifies significantly co-occurring potential latent drivers. We observe 155 same gene double mutations of which 140 individual components are cataloged as latent drivers. Evaluation of cell lines and patient-derived xenograft response data to drug treatment indicate that in certain genes double mutations may have a prominent role in increasing oncogenic activity, hence obtaining a better drug response, as in PIK3CA. Taken together, our comprehensive analyses indicate that same-gene double mutations are exceedingly rare phenomena but are a signature for some cancer types, e.g., breast, and lung cancers. The relative rarity of doublets can be explained by the likelihood of strong signals resulting in oncogene-induced senescence, and by doublets consisting of non-identical single residue components populating the background mutational load, thus not identified.

Plasma proteomics identify potential severity biomarkers from COVID-19 associated network.

PURPOSE: Coronavirus disease 2019 (COVID-19) continues to threaten public health globally. Severe acute respiratory coronavirus type 2 (SARS-CoV-2) infection-dependent alterations in the host cell signaling network may unveil potential target proteins and pathways for therapeutic strategies. In this study, we aim to define early severity biomarkers and monitor altered pathways in the course of SARS-CoV-2 infection. EXPERIMENTAL DESIGN: We systematically analyzed plasma proteomes of COVID-19 patients from Turkey by using mass spectrometry. Different severity grades (moderate, severe, and critical) and periods of disease (early, inflammatory, and recovery) are monitored. Significant alterations in protein expressions are used to reconstruct the COVID-19 associated network that was further extended to connect viral and host proteins. RESULTS: Across all COVID-19 patients, 111 differentially expressed proteins were found, of which 28 proteins were unique to our study mainly enriching in immunoglobulin production. By monitoring different severity grades and periods of disease, CLEC3B, MST1, and ITIH2 were identified as potential early predictors of COVID-19 severity. Most importantly, we extended the COVID-19 associated network with viral proteins and showed the connectedness of viral proteins with human proteins. The most connected viral protein ORF8, which has a role in immune evasion, targets many host proteins tightly connected to the deregulated human plasma proteins. CONCLUSIONS AND CLINICAL RELEVANCE: Plasma proteomes from critical patients are intrinsically clustered in a distinct group than severe and moderate patients. Importantly, we did not recover any grouping based on the infection period, suggesting their distinct proteome even in the recovery phase. The new potential early severity markers can be further studied for their value in the clinics to monitor COVID-19 prognosis. Beyond the list of plasma proteins, our disease-associated network unravels altered pathways, and the possible therapeutic targets in SARS-CoV-2 infection by connecting human and viral proteins. Follow-up studies on the disease associated network that we propose here will be useful to determine molecular details of viral perturbation and to address how the infection affects human physiology.

Structural coverage of the human interactome.

Complex biological processes in cells are embedded in the interactome, representing the complete set of protein-protein interactions. Mapping and analyzing the protein structures are essential to fully comprehending these processes' molecular details. Therefore, knowing the structural coverage of the interactome is important to show the current limitations. Structural modeling of protein-protein interactions requires accurate protein structures. In this study, we mapped all experimental structures to the reference human proteome. Later, we found the enrichment in structural coverage when complementary methods such as homology modeling and deep learning (AlphaFold) were included. We then collected the interactions from the literature and databases to form the reference human interactome, resulting in 117 897 non-redundant interactions. When we analyzed the structural coverage of the interactome, we found that the number of experimentally determined protein complex structures is scarce, corresponding to 3.95% of all binary interactions. We also analyzed known and modeled structures to potentially construct the structural interactome with a docking method. Our analysis showed that 12.97% of the interactions from HuRI and 73.62% and 32.94% from the filtered versions of STRING and HIPPIE could potentially be modeled with high structural coverage or accuracy, respectively. Overall, this paper provides an overview of the current state of structural coverage of the human proteome and interactome.

Comparative biological network analysis for differentially expressed proteins as a function of bacilysin biosynthesis in Bacillus subtilis.

The Gram-positive bacterium Bacillus subtilis produces a diverse range of secondary metabolites with different structures and activities. Among them, bacilysin is an enzymatically synthesized dipeptide that consists of L-alanine and L-anticapsin. Previous research by our group has suggested bacilysin's role as a pleiotropic molecule in its producer, B. subtilis PY79. However, the nature of protein interactions in the absence of bacilysin has not been defined. In the present work, we constructed a protein-protein interaction subnetwork by using Omics Integrator based on our recent comparative proteomics data obtained from a bacilysin-silenced strain, OGU1. Functional enrichment analyses on the resulting networks pointed to certain putatively perturbed pathways such as citrate cycle, quorum sensing and secondary metabolite biosynthesis. Various molecules, which were absent from the experimental data, were included in the final network. We believe that this study can guide further experiments in the identification and confirmation of protein-protein interactions in B. subtilis.

Functional stratification of cancer drugs through integrated network similarity.

Drugs not only perturb their immediate protein targets but also modulate multiple signaling pathways. In this study, we explored networks modulated by several drugs across multiple cancer cell lines by integrating their targets with transcriptomic and phosphoproteomic data. As a result, we obtained 236 reconstructed networks covering five cell lines and 70 drugs. A rigorous topological and pathway analysis showed that chemically and functionally different drugs may modulate overlapping networks. Additionally, we revealed a set of tumor-specific hidden pathways with the help of drug network models that are not detectable from the initial data. The difference in the target selectivity of the drugs leads to disjoint networks despite sharing a similar mechanism of action, e.g., HDAC inhibitors. We also used the reconstructed network models to study potential drug combinations based on the topological separation and found literature evidence for a set of drug pairs. Overall, network-level exploration of drug-modulated pathways and their deep comparison may potentially help optimize treatment strategies and suggest new drug combinations.

Context dependent isoform specific PI3K inhibition confers drug resistance in hepatocellular carcinoma cells.

BACKGROUND: Targeted therapies for Primary liver cancer (HCC) is limited to the multi-kinase inhibitors, and not fully effective due to the resistance to these agents because of the heterogeneous molecular nature of HCC developed during chronic liver disease stages and cirrhosis. Although combinatorial therapy can increase the efficiency of targeted therapies through synergistic activities, isoform specific effects of the inhibitors are usually ignored. This study concentrated on PI3K/Akt/mTOR pathway and the differential combinatory bioactivities of isoform specific PI3K-α inhibitor (PIK-75) or PI3K-ß inhibitor (TGX-221) with Sorafenib dependent on PTEN context. METHODS: The bioactivities of inhibitors on PTEN adequate Huh7 and deficient Mahlavu cells were investigated with real time cell growth, cell cycle and cell migration assays. Differentially expressed genes from RNA-Seq were identified by edgeR tool. Systems level network analysis of treatment specific pathways were performed with Prize Collecting Steiner Tree (PCST) on human interactome and enriched networks were visualized with Cytoscape platform. RESULTS: Our data from combinatory treatment of Sorafenib and PIK-75 and TGX-221 showed opposite effects; while PIK-75 displays synergistic effects on Huh7 cells leading to apoptotic cell death, Sorafenib with TGX-221 display antagonistic effects and significantly promotes cell growth in PTEN deficient Mahlavu cells. Signaling pathways were reconstructed and analyzed in-depth from RNA-Seq data to understand mechanism of differential synergistic or antagonistic effects of PI3K-α (PIK-75) and PI3K-ß (TGX-221) inhibitors with Sorafenib. PCST allowed as to identify AOX1 and AGER as targets in PI3K/Akt/mTOR pathway for this combinatory effect. The siRNA knockdown of AOX1 and AGER significantly reduced cell proliferation in HCC cells. CONCLUSIONS: Simultaneously constructed and analyzed differentially expressed cellular networks presented in this study, revealed distinct consequences of isoform specific PI3K inhibition in PTEN adequate and deficient liver cancer cells. We demonstrated the importance of context dependent and isoform specific PI3K/Akt/mTOR signaling inhibition in drug resistance during combination therapies. ( https://github.com/cansyl/Isoform-spesific-PI3K-inhibitor-analysis ).

Artificial intelligence based methods for hot spot prediction.

Proteins interact through their interfaces to fulfill essential functions in the cell. They bind to their partners in a highly specific manner and form complexes that have a profound effect on understanding the biological pathways they are involved in. Any abnormal interactions may cause diseases. Therefore, the identification of small molecules which modulate protein interactions through their interfaces has high therapeutic potential. However, discovering such molecules is challenging. Most protein-protein binding affinity is attributed to a small set of amino acids found in protein interfaces known as hot spots. Recent studies demonstrate that drug-like small molecules specifically may bind to hot spots. Therefore, hot spot prediction is crucial. As experimental data accumulates, artificial intelligence begins to be used for computational hot spot prediction. First, we review machine learning and deep learning for computational hot spot prediction and then explain the significance of hot spots toward drug design.

Mechanism of activation and the rewired network: New drug design concepts.

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Computational approaches leveraging integrated connections of multi-omic data toward clinical applications.

In line with the advances in high-throughput technologies, multiple omic datasets have accumulated to study biological systems and diseases coherently. No single omics data type is capable of fully representing cellular activity. The complexity of the biological processes arises from the interactions between omic entities such as genes, proteins, and metabolites. Therefore, multi-omic data integration is crucial but challenging. The impact of the molecular alterations in multi-omic data is not local in the neighborhood of the altered gene or protein; rather, the impact diffuses in the network and changes the functionality of multiple signaling pathways and regulation of the gene expression. Additionally, multi-omic data is high-dimensional and has background noise. Several integrative approaches have been developed to accurately interpret the multi-omic datasets, including machine learning, network-based methods, and their combination. In this review, we overview the most recent integrative approaches and tools with a focus on network-based methods. We then discuss these approaches according to their specific applications, from disease-network and biomarker identification to patient stratification, drug discovery, and repurposing.

Ischemic Heart Disease Selectively Modifies the Right Atrial Appendage Transcriptome.

Background: Although many pathological changes have been associated with ischemic heart disease (IHD), molecular-level alterations specific to the ischemic myocardium and their potential to reflect disease severity or therapeutic outcome remain unclear. Currently, diagnosis occurs relatively late and evaluating disease severity is largely based on clinical symptoms, various imaging modalities, or the determination of risk factors. This study aims to identify IHD-associated signature RNAs from the atrial myocardium and evaluate their ability to reflect disease severity or cardiac surgery outcomes. Methods and Results: We collected right atrial appendage (RAA) biopsies from 40 patients with invasive coronary angiography (ICA)-positive IHD undergoing coronary artery bypass surgery and from 8 patients ICA-negative for IHD (non-IHD) undergoing valvular surgery. Following RNA sequencing, RAA transcriptomes were analyzed against 429 donors from the GTEx project without cardiac disease. The IHD transcriptome was characterized by repressed RNA expression in pathways for cell-cell contacts and mitochondrial dysfunction. Increased expressions of the CSRNP3, FUT10, SHD, NAV2-AS4, and hsa-mir-181 genes resulted in significance with the complexity of coronary artery obstructions or correlated with a functional cardiac benefit from bypass surgery. Conclusions: Our results provide an atrial myocardium-focused insight into IHD signature RNAs. The specific gene expression changes characterized here, pave the way for future disease mechanism-based identification of biomarkers for early detection and treatment of IHD.

Epithelial Wnt secretion drives the progression of inflammation-induced colon carcinoma in murine model.

Colon cancer is initiated by stem cells that escape the strict control. This process is often driven through aberrant activation of Wnt signaling by mutations in components acting downstream of the receptor complex that unfetter tumor cells from the need for Wnts. Here we describe a class of colon cancer that does not depend on mutated core components of the Wnt pathway. Genetically blocking Wnt secretion from epithelial cells of such tumors results in apoptosis, reduced expression of colon cancer markers, followed by enhanced tumor differentiation. In contrast to the normal colonic epithelium, such tumor cells autosecrete Wnts to maintain their uncontrolled proliferative behavior. In humans, we determined certain cases of colon cancers in which the Wnt pathway is hyperactive, but not through mutations in its core components. Our findings illuminate the path in therapy to find further subtypes of Wnt-dependent colon cancer that might be responsive to Wnt secretion inhibitors.

Performance Assessment of the Network Reconstruction Approaches on Various Interactomes.

Beyond the list of molecules, there is a necessity to collectively consider multiple sets of omic data and to reconstruct the connections between the molecules. Especially, pathway reconstruction is crucial to understanding disease biology because abnormal cellular signaling may be pathological. The main challenge is how to integrate the data together in an accurate way. In this study, we aim to comparatively analyze the performance of a set of network reconstruction algorithms on multiple reference interactomes. We first explored several human protein interactomes, including PathwayCommons, OmniPath, HIPPIE, iRefWeb, STRING, and ConsensusPathDB. The comparison is based on the coverage of each interactome in terms of cancer driver proteins, structural information of protein interactions, and the bias toward well-studied proteins. We next used these interactomes to evaluate the performance of network reconstruction algorithms including all-pair shortest path, heat diffusion with flux, personalized PageRank with flux, and prize-collecting Steiner forest (PCSF) approaches. Each approach has its own merits and weaknesses. Among them, PCSF had the most balanced performance in terms of precision and recall scores when 28 pathways from NetPath were reconstructed using the listed algorithms. Additionally, the reference interactome affects the performance of the network reconstruction approaches. The coverage and disease- or tissue-specificity of each interactome may vary, which may result in differences in the reconstructed networks.

Epitranscriptomics of Ischemic Heart Disease-The IHD-EPITRAN Study Design and Objectives.

Epitranscriptomic modifications in RNA can dramatically alter the way our genetic code is deciphered. Cells utilize these modifications not only to maintain physiological processes, but also to respond to extracellular cues and various stressors. Most often, adenosine residues in RNA are targeted, and result in modifications including methylation and deamination. Such modified residues as N-6-methyl-adenosine (m6A) and inosine, respectively, have been associated with cardiovascular diseases, and contribute to disease pathologies. The Ischemic Heart Disease Epitranscriptomics and Biomarkers (IHD-EPITRAN) study aims to provide a more comprehensive understanding to their nature and role in cardiovascular pathology. The study hypothesis is that pathological features of IHD are mirrored in the blood epitranscriptome. The IHD-EPITRAN study focuses on m6A and A-to-I modifications of RNA. Patients are recruited from four cohorts: (I) patients with IHD and myocardial infarction undergoing urgent revascularization; (II) patients with stable IHD undergoing coronary artery bypass grafting; (III) controls without coronary obstructions undergoing valve replacement due to aortic stenosis and (IV) controls with healthy coronaries verified by computed tomography. The abundance and distribution of m6A and A-to-I modifications in blood RNA are charted by quantitative and qualitative methods. Selected other modified nucleosides as well as IHD candidate protein and metabolic biomarkers are measured for reference. The results of the IHD-EPITRAN study can be expected to enable identification of epitranscriptomic IHD biomarker candidates and potential drug targets.

Structural analysis of mammalian protein phosphorylation at a proteome level.

Phosphorylation is an essential post-translational modification for almost all cellular processes. Several global phosphoproteomics analyses have revealed phosphorylation profiles under different conditions. Beyond identification of phospho-sites, protein structures add another layer of information about their functionality. In this study, we systematically characterize phospho-sites based on their 3D locations in the protein and establish a location map for phospho-sites. More than 250,000 phospho-sites have been analyzed, of which 8,686 sites match at least one structure and are stratified based on their respective 3D positions. Core phospho-sites possess two distinct groups based on their dynamicity. Dynamic core phosphorylations are significantly more functional compared with static ones. The dynamic core and the interface phospho-sites are the most functional among all 3D phosphorylation groups. Our analysis provides global characterization and stratification of phospho-sites from a structural perspective that can be utilized for predicting functional relevance and filtering out false positives in phosphoproteomic studies.

Quantitative Proteomics Identifies Secreted Diagnostic Biomarkers as well as Tumor-Dependent Prognostic Targets for Clear Cell Renal Cell Carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the third most common and most malignant urological cancer, with a 5-year survival rate of 10% for patients with advanced tumors. Here, we identified 10,160 unique proteins by in-depth quantitative proteomics, of which 955 proteins were significantly regulated between tumor and normal adjacent tissues. We verified four putatively secreted biomarker candidates, namely, PLOD2, FERMT3, SPARC, and SIRPα, as highly expressed proteins that are not affected by intratumor and intertumor heterogeneity. Moreover, SPARC displayed a significant increase in urine samples of patients with ccRCC, making it a promising marker for the detection of the disease in body fluids. Furthermore, based on molecular expression profiles, we propose a biomarker panel for the robust classification of ccRCC tumors into two main clusters, which significantly differed in patient outcome with an almost three times higher risk of death for cluster 1 tumors compared with cluster 2 tumors. Moreover, among the most significant clustering proteins, 13 were targets of repurposed inhibitory FDA-approved drugs. Our rigorous proteomics approach identified promising diagnostic and tumor-discriminative biomarker candidates which can serve as therapeutic targets for the treatment of ccRCC. IMPLICATIONS: Our in-depth quantitative proteomics analysis of ccRCC tissues identifies the putatively secreted protein SPARC as a promising urine biomarker and reveals two molecular tumor phenotypes.